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We theoretically investigate negative differential resistance (NDR) in graphene double barrier 
resonant tunneling diodes within the Landauer-Biittiker formalism by considering a realistic linear 
voltage drop. The calculated current-voltage characteristic indicates that, due to the competition 
between hole-to-electron transport, Klein paradox, and resonant tunneling, the NDR is substantially 
suppressed and only presents for appropriate structural parameters. Remarkably, the first NDR 
operation window is always pinned by a robust kink effect to around the Fermi energy, which 
may open an interesting window for designing NDR based graphene devices in a controllable way. 
Moreover, the peak-to-valley current ratio can be enhanced epitaxially and parabolically by opening 
a band gap and adjusting structural asymmetries, respectively. 



I. INTRODUCTION 

The negative differential resistance (NDR) is a fun- 
damental phenomenon in physics and can be utilized to 
generate bistable switching circuits and electronic oscilla- 
tors It has been observed in various systems, ranging 
from gaseous medium [21 via chalcogenide glasses [3 to 
organic semiconductors j and conductive polymers 
Recently, significant efforts have been devoted to 

the study of NDR in graphene, a monolayer of sp^ bonded 
carbon atoms that has attracted much attentions since its 
discovery [l^]. Besides its occurrence in graphene single 
barrier diodes (SBDs) the NDR has also been the- 

oretically reported in zigzag graphene nanoribbons 0, 
armchair graphene nanoribbon s fl0|. p TI. and numerus 



graphene nanoribbon junctions Tj governed by the 
mechanisms of chiral selective tunneling, electrodes dop- 
ing or strain, and suppression of the coherent transition, 
respectively. Very recently, it is suggested that the NDR 
can be obtained in armchair graphene superlattices (iBj 
by band gaps, miniband conductance, and Wannier-Stark 
ladder and rather surprisingly in graphene three terminal 
field-effect transistors [l^ by just considering classical 
ambipolar transport. 

It is well-known that Fabry-Perot-type interference is 
an alternatively fundamental mechanism for NDR and 
plays a dominant role on the NDR in common semi- 
conductor based resonant tunneling diodes (RTDs) fisj . 
Surprisingly, this mechanism has not been explored in 
graphene so far except a rather recent study on semi- 
conducting armchair graphene nanoribbon superlattices 
[Tsj in which, however, other mechanisms as mentioned 
above also contribute to the NDR features thus signif- 
icantly obscure the resonant tunneling mechanism. On 
the other hand, in those above works based on analytic 
models [l-d, [13 the voltage drop produced by the bias 
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has been approximated by a steplike one of the average 
value. (See, Figs. 1, 3, and 1 in Refs. Q, 0, and [l5l |. 
respectively.) This approximation greatly simplifies the 
calculations (see below), but actually describes a bias- 
dependent nonequlibrium barrier which may draw a false 
conclusion of the NDR feature. 

In present work, we for the first time consider a realistic 
linear voltage drop and investigate the NDR in graphene 
double barrier (DB) RTDs. The said DB RTDs can be 
fabricated by patterning two top gates [l^ on a short 
and sufficient wide graphene strip [20|, thus the edge ef- 
fect can be ignored and the NDR feature solely relies on 
the resonant tunneling mechanism. To investigate the I- 
V characteristic, we analytically solve the envelope wave 
functions in a rotated pseudospin space. We find that the 
NDR in graphene DB RTDs has been substantially sup- 
pressed comparing with that in common semiconductors, 
and it only appears with appropriate structural parame- 
ters. We attribute such a suppression to hole-to-electron 
transport and the Klein paradox, and propose a way to 
enhance the strength of NDR (defined as the peak-to- 
valley current ratio (PVR)) by introducing a band gap in 
the linear dispersion of graphene or adjusting the struc- 
tural parameters. Remarkably, our calculations show 
that the NDR is pinned by a so-called kink effect and 
the first operation window always appears around the 
Fermi energy [2l] | . This phenomenon has not been found 
in other type NDR devices and may provide an interest- 
ing opportunity for designing and utilizing NDR based 
devices in a controllable way. Moreover, our results also 
indicate that the NDR is absent in graphene SBDs as 
long as the linear voltage drop is taken into account. 

The paper is organized as follows. In Sec. II we intro- 
duce the theoretical model and formula that we use in 
order to calculate the I-V curves. In Sec. Ill we present 
the TV and DC characteristics for graphene SBDs con- 
cerning the linear voltage drop. In Sec. IV we exhibit 
the kink-effect pinning NDR in graphene DB RTDs and 
in Sec. V we investigate the role of band gaps and struc- 
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FIG. 1. (color online) Model construction for I-V characteris- 
tic of a graphene DB RTD, under (a) no bias and (b) a finite 
bias Vb- The diode contains a buffer (width wi,), two barriers 
(width uii(2) and height Vti(2)), and a well (width d) regions. 
For (b) a linear voltage drop appears and the transport can 
be divided into three regimes (I, II, and III) defined by Ep 
and Vb- (c) The rotation we make in the pseudospin space. 



tural asymmetries on the NDR strength. A summary will 
be given in Sec. VI. 



II. MODEL AND FORMULA 

Fig. 1(a) shows a schematic of a graphene DB RTD 
under zero bias. To investigate its intrinsic I-V charac- 
teristic, we restrict our study to the ideal contacts case 
which in practice can be achieved by usinga high-K (di- 
electric constant) material as electrodes [22|. The con- 
tact resistance and the doping effect by electrodes then 
can be ignored. The Fermi energy Ep (respective to 
the Dirac point Ejj = 0) and the carrier concentration 
(oc \Ef\) [i31 in the stripe are tuned by a back gate volt- 
age. The potential barriers (Vti(2)) can be formed by 
top gates [13. We adopt rectangular potential barri- 
ers, since smooth barriers considering the electric field 
on the interfaces can be regarded as rectangular ones 
with the same widths but effective heights determined 
by the smoothness As sketched in Fig. 1(b), when 
a bias voltage (VJ,) is applied, a linear voltage drop will 
be formed between the two electrodes due to a uniform 
electric field and a nest current will be produced by the 
holes or electrons with energies E S [Ep — eVb, Ep]. The 
charge transport can be classified into three regimes de- 
pending on the values of Ep and Vb as denoted in Fig. 



1(b). In regime I (0 < ii^ < Ep), it is electron-to-electron 
transport; while in both regimes II (— eVb/2 < E < 0) 
and III {Ep — eVb < E < — eVt/2), it is hole-to-electron 
transport. The difference between the latter two regimes 
is that in regime III only holes from the left electrode 
within critical incident angles =Fsin"^(l + eVb/E) can 
transport to electrons in the right electrode. One can 
see that when eVb < Ep, only regime I is possible; when 
Ep < eVb < 2Ep, regime II is also present; while for 
eVb > 2Ep all these three regimes make contributions to 
the nest current. 

At OK the nest current under a bias Vb can be given 
by the Landauer-Biittiker formalism [23] 

I{Vb) = hr I ' T{e,e,Vb)\t\cosedede, (1) 

where Iq = ^evpWy/ is a current unit. The 
factor 4 comes from the spin and valley degeneracies, 
vp w 10^to/s is the Fermi velocity, and Wy is the width 
of the graphene strip, is a characteristic length of 
the structure which is adopted as 40 nm to ensure the 
electron density of states coinciding with a true system 
and a coherent transport regime at room temperature. 
Ep and Vb in the integral are written in the dimen- 
sionlcss form; sec the next paragraph for details. To 
obtain the value of the current, one need to calculate 
the transmission probabilities (T) at i? G [Ep — eVb, Ep] 
first. This can be obtained by matching spinor wavcfunc- 
tions {^{x,y)) in different regions at potential bound- 
aries to ensure the conservation of the local current den- 
sity (j(r) = Aevp'^\v)a'^{Y)). However, direct deriva- 
tions from the Dirac equation with a horizontal electric 
field result in a unsolvable two-order differential equa- 
tion. In Refs. [Gt-lsl-fTst this dilemma has be circumvented 
by approximating the linear voltage drop with an aver- 
age value. However, such an approximation describes a 
nonequlibrium barrier which may induce a nonphysical 
current. 

Here, we perform a rotation of the Dirac equation by 
7r/2 around the axis in the pseudospin space (see. Fig. 
1(c)). After rotation the Hamiltonian becomes 

H = vp[a^{p.^+eA^)+ay{py+eAy)-a^/S\+eV{x)l, (2) 

and the two components of the spinor ('0^) relate to the 
ones before rotation (ijA-^) by -0^ = (iV'^ + '/'^)/'\/2- 
Here il/''/-^{x, y) is the envelope function at sublattice site 
A/B of the graphene sheet, cr = {<Jx,o'y,crz) are Pauli's 
matrices, I is the 2x2 identity matrix, p = {px,Py)'^ = 
{hkx, hky)^ is the momentum operator, A — {A^, Ay)'^ is 
a possible magnetic vector potential, A is a tunable band 
gap up to several hundreds of meV achieved through a 
controllable doping [251, and V{x) is the position de- 
pendent electrostatic potential. In terms of do all these 
parameters can be expressed in a dimensionless form: 
X = x/do, k = kdo, e = E/Eq, S = A/Eq, and v{x) = 
eV{x)/EQ with an energy unit Eq = hvp/do « 16.44 
mcV. 
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Prom Eq. (2) we get coupled one-order differential 
equations [id / dx±{e — v{x))]il:^ = {—d/dy^S)^^ . Tak- 
ing into account solutions of ip^{x,y) — /^(x)e*'^!'^ due 
to the translation invariance in the y-axis, we get a de- 
coupled two-order differential equation 

± *[^' K^)] + - '^C^)]' f^i^) = 0' (3) 

where [d/dx,v{x)] is a commutator due to the noncom- 
mutativity of Px and x [1^ and = ky + S"^. In the 
electrode regions iv{x) =const), this equation recovers 
the one before performing rotation, [cP/dx"^ + ^ 'f^)^ ^ 
q^]f^{x) = 0. The direct wave function solutions from 
this equation should be composed as mentioned above. In 
the buffer and well (barrier) regions, v{x) = —ax{+vt), 
where a = eVb/L and L is the total width between the 
electrodes. Now, the wave functions of Eq. (3) can be 
solved as 

(4) 

where F = D[-l + iq^/2a,{l + ax{-Vt)) / y/a]^ 

G = {l + i)y/aq-^ D[iq^ /2a, {l + i){e + ax{-vt))l and 
the superscript stands for complex conjugate. Here D is 
the Weber parabolic cylinder function. F and G have 
properties of a right-going wave function as the amplitude 
of the classical trajectory at the right turning point (the 
points where kx{x) = •\/(e — v{x)Y — 9^ becomes zero) 
is about e"'^' '° times of that at the left turning point. 
Since we have got wave functions in the electrodes and 
intermediate regions, we can obtain the transmission co- 
efficient {t) by the standard transfer-matrix method [l^] . 
The transmission probability reads T — ^^^f^jjj^^^jp-\t\'^ 
for fc^^ > and fc^^^ > 0, and T — Q otherwise, where 
kxL(R) is the value of kx at the left (right) electrode. 
Then the current can be calculated by Eq. (1). 

III. I-V CHARACTERISTIC IN GRAPHENE 
SBDS 

To evaluate the results taking into account a realistic 
linear voltage drop, and more importantly, to obtain a 
basic picture for biased transport in graphene in absence 
of resonant tunneling, we'd like to revisit the NDR in 
graphene SBDs first. Fig. 2 shows the I-V and differen- 
tial conductance (DC) curves for a graphene SBD with 
Wb = 1, w = 1, and 7;t = 2 at various back gate volt- 
ages. As is seen, for zero back gate (i.e., Ep = 0), the 
I-V curve is superlinear (/ oc Vjf ) for VJ, < Vt and nearly 
linear (/ cx Vf,) for Vb > Vt- This is a result of that the 
current is totally due to the hole-to-electron Zener-Klein 
tunneling. As long as the Fermi energy is lifted, the DC 
displays totally different behaviors in the three transport 
regimes (see. Fig. 2(b)). Take cf = 1 (black dotted line) 
as an example. In regime I, the DC decreases with in- 
creasing Vb and reaches Gm exactly at eVb = Ep. This is 
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FIG. 2. (color online) (a) I-V and (b) DC characteristics for 
a graphene SBD with parameters of uij, = 1, w = 1, vt = 2. 
Along the arrow, eF=0, 1, 1.5, 2, 2.5, and 3, respectively. In 
(b) the DC is also shown for eF=0.25 (line with dots) and 0.5 
(line with circles) in small ranges. The conductance unit is 
Go = 4e'^Wy/2nhdo. 



because the bias-induced electron concentration (approx- 
imately oc E~eVb) which contributes to the DC decreases 
with the bias. Further increasing eVb to 2Ep (regime II), 
the DC increases with increasing Vb and reaches a maxi- 
mum at Vh = 2Ep. In this regime, the bias-induced hole 
concentration (approximately oc eVb — E) increases with 
bias. When eVb exceeds 2Ep (regime III), increasing bias 
will increase the hole concentration but decrease the DC 
as only holes with incident angles smaller than the criti- 
cal value contribute to the transport. As a result of such 
an ambipolar transport by both electron and hole, the 
I-V curve shows a kink effect exactly at eVh = Ep, which 
has been observed by recent experiments [28| . 

Interestingly, Gm increases with increasing Ep (the 
ep = 2 case is an exception maybe due to it equals exactly 
to the barrier potential) and is positive even for a rather 
small Ep. In other words, there is no NDR effect in 
graphene-based SBDs. This conclusion is in strong con- 
trast to the one obtained by Dragomans @ and Nam Do 
et al. 0, Q . The artificial current induced by the steplike 
voltage drop in their works may be responsible for the an- 
nounced NDR effect. Note, the reduction in T may give 
rise to NDR in SBDs based on common semiconductors 
(such as Hgi_:,Cd^Tc/CdTc ^ and AlGaSb/InAs ^). 
Here the situation substantially changes for graphene. 
Due to the transport regimes II and HI which are absent 
in common semiconductor cases, the number of tunneling 
carriers (both electrons and holes) increase rather than 
saturate with increasing bias. Then the reduction in T is 
obscured and the NDR effect is absent in graphene SBDs. 
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FIG. 3. (color online) (a) I-V and (b) DC characteristics for 
a graphene DB RTD with parameters of Wh = 1, ti'i(2) = 1, 
d = 5, and Vt\(2) = 1- Along the arrow, e_F=0.7, 0.8, 0.9, 
1.0, and 1.1. Insert in (a): I-V characteristic for a 2DEG 
DB RTD with the same parameters (note now Eq = ft^/2mdo 
with m the electron mass in 2DEG). The thin line shows I-V 
characteristic for a 2DEG SBD with Ef = 0.2Eo. Insert in 
(b): the minimum DC as a function of the Fermi energy. 



IV. NDR IN GRAPHENE DB RTDS 

We now seek NDR in graphene DB RTDs and investi- 
gate its correlation with the resonant tunneling. Fig. 3 
shows the I-V and DC characteristics for a graphene DB 
RTD with Wb = 1, wi = W2 = d = 5, and vn = f t2 = 1 
at various back gating. While exhibiting similar trend as 
the ones for a graphene SBD, the I-V curves here for 
a graphene DB RTD display obvious ripples and corre- 
spondingly the DC oscillates with the bias. This is a 
result of the alternate enhancement and suppression of 
T respectively around and between resonant tunnelings 
(see, Fig. 4(a)). For the latter case, the DC reaches 
local minimums Gim at some biases. One can see from 
Fig. 3(b) that at biases sufficiently low or high compared 
with the Fermi energies, Gim moves to a lower bias for 
a higher Fermi energy (see, the groups around Vb =0.2, 
1.7, and 2.8; the one around Vb = 1.2 for = 0.7 is 
excluded in this analysis). However, the lowest Gim (i-e., 
Gm) always locates exactly at eVb = Ep- It first de- 
creases with increasing Ep until reaching the minimum 
DC of the structure, and then increases with increasing 
Ep (see, insert in Fig. 3(b)). Note, the minimum DC of 
the structure depends only on the structural parameters 
and for the considered structure, the minimum DC ap- 
pears at Ep = 0.9. In total, the G„i is always positive, 
which means that there is no NDR in such a graphene 
DB RTD. 

This is an interesting conclusion since the NDR always 
occurs for common semiconductor based DB RTDs (see 
inset in Fig. 3(a)). It is well known that in the semicon- 
ductor RTDs, NDR occurs when the suppression region 
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FIG. 4. (color online) Difference of transmission probabilities 
between a unbiased DB RTD and constituted single barrier 
for (a) graphene and (b) 2DEG with the same parameters of 
wi(2) = 1, d = 5, and Vti(2) = 1- 



of T between two adjacent resonant tunneling peaks (i.e., 
blue regions in Fig. 4(b)) enters the integration window 
[Ep — eVb,Ep] of the current [l^- Meanwhile, one can 
see from Eq. (1) that the contribution of the suppression 
regions to the NDR mainly comes from ones near = 
due to the cos 9 factor. However, in graphene these re- 
gions are significantly reduced. Due to the Klein paradox, 
the suppression region between quasibound states (equiv- 
alcntly, resonant tunnelings) cannot form especially for 
small incident angles (see. Fig. 4(a)). Moreover, as indi- 
cated in the above graphene SBD calculations, when eVb 
exceeds Ep the hole-to-electron tunneling contributes a 
positive DC (which increases with increasing bias) hence 
suppresses the NDR feature in graphene. Therefore, the 
NDR effect in graphene DB RTDs is a competition be- 
tween hole-to-electron transport, Klein paradox, and res- 
onant tunneling. 

Our calculations further indicate that, the absence of 
NDR can be overcome by enhancing the resonant tunnel- 
ing with more appropriate structural parameters. Specif- 
ically, the less the number of quasibound states (which 
approximately equals to the value of Vtd), the stronger 
the contribution of the resonant tunneling. Fig. 5 shows 
the results for a structure of Wb = 1/2, wi = W2 = 1/2, 
d = 1/2, and vn = Vt2 = 2. One can see that Gm be- 
comes negative for ep ^ 2.6 — 3.2. Therefore the NDR 
can be obtained in graphene DB RTDs with proper struc- 
tural parameters. It is noted in Fig. 5(b) that, even in 
the presence of NDR feature. Gm is also locked exactly at 
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FIG. 5. (color online) (a) I-V and (b) DC characteristics for 
a graphene DBD with parameters of wt = 1/2, Wn2) = 1/2, 
d = 1/2, and Vti(2) = 2. Along the arrow, 6^=2.4, 2.6, 2.8, 
3.0, and 3.2. 



eVb = Ep, implying that the kink effect is robust against 
the resonant tunneling. This means that the first NDR 
operation window can be chosen or tuned easily by a 
single parameter, i.e., the Fermi energy. 



V. EFFECT OF BAND GAPS AND 
STRUCTURAL ASYMMETRIES 

However, the value of PVR is rather small (about 1.02 
for the strongest case ep = 2.8) for these symmetry struc- 
tures based on gaplcss graphene. Fig. 6(a) shows the 
I-V characteristic for a graphene DB RTD with differ- 
ent band gaps. As is seen, when there is an induced 
band gap in graphene (which should be restricted to 
not exceed the value of Fermi energy), the current de- 
creases for any given bias. This is a result of the re- 
duction of number of electron/holes from the left elec- 
trode (i.e., the ones with |ecos0| < 6 cannot propagate 
freely) and the suppression of T by the band gap. When 
the band gap is sufficiently large (say S > 1.5), the I-V 
characteristic at low bias region even becomes superlin- 
ear and a second NDR operation window with a much 
weaker strength appears at higher bias (see. Figs. 7 and 
8). However, the PVR for the first operation window 
does increase with increase S (see. Fig. 6(b)). The in- 
crease follows an exponential function of the band gap 
(PVR sa 0.9832 + 0M13c^-^^^ for the considered case). 
The underlying physics is that, the presence of band gap 
increases the modulus of the image longitudinal wave vec- 
tor in the barriers {kj:{x) = [6^ + ky — {e — v{x))'^]~^^^), 
which makes a stronger suppression of T between adja- 
cent resonant tunnelings hence results in a bigger PVR. 

We now study the effect of the structural asymme- 
tries of DB RTDs, i.e., the two potential barriers have 
different widths and/or heights. To give a quantitative 



FIG. 6. (color online) (a) I-V characteristics for the graphene 
DB RTD considered in Fig. 5 with different band gaps. Along 
the arrow, A/Eo=0, 0.5, 1.0, 1.5, 2.0, and 2.5 (enlarged by 5 
times for clearness), (b) The peak-to- valley current ratio and 
(c) the biases for the current peak, kink, and valley points as 
a function of the induced gap. 



evaluation, here we suppose that the height of the left 
(right) barrier is lower (higher) than a referenced sym- 
metric one by a ratio of Xh = {vt~vti)/vt = {vt2~vt)/vt, 
which can be defined as an asymmetry degree in bar- 
rier heights (ADBH). Fig. 7(a) shows I-V characteristics 
for graphene DB RTDs with different ADBHs and cor- 
responding PVR's and biases for the peaks, kinks, and 
valleys are plotted in Fig. 7(b) and 7(c), respectively. 
As is seen, the I-V curves exhibit regular variations with 
increasing ADBH, not only in gradually decreasing bias 
positions for the peak, kink, and valley but also in a 
parabolic-like ADBH dependent PVR. The maximum 
PVR is achieved at about Xh = 0.2 (i.e., Vti:vt2='2:3). 
Note, Xh = stands for a symmetric DB RTD, for which 
the resonant tunneling does not achieve the optimal one 
under bias [3ll | ; while A/i = 1 means actually a single bar- 
rier diode, for which the NDR never happen as we have 
shown in Section II. This qualitatively explains why the 
maximum PVR happens at an intermediate ADBH be- 
tween and 1. 

Fig. 8(a) shows I-V curves for different asymmetry 
degrees in barrier widths (ADBWs), which are similarly 
defined as A„, = {w — wi)/w = {w2 — w)/w. It is seen 
that, the low bias behaviors almost coincide no matter 
how big the ADBW is. Different from the changes due 
to asymmetry in barrier heights, where all the biases 
for the peak, kink, and valley decrease with increasing 
ADBH. here the biases for the peak and kink increase 
slightly with increasing ADBW, while the one for the val- 
ley shows a weak parabolic-like behavior (see. Fig. 8(c)). 
More importantly, the PVR also shows a parabolic-like 
dependence on the ADBW with the maximum also at 
about A„ = 0.2 {wi:w2=2:S) (see. Fig. 8(b)). This sim- 
ilarity in the PVR behaviors implies the similar effects 
of asymmetries in barrier heights and widths. Further 
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parameters. We have also proposed a way to enhance 
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FIG. 7. (color online) (a) I-V characteristics for the graphene 
DB RTD considered in Fig. 5 with 5 = 2.5 and different 
ADBHs. Along the arrow, Xh=0, 0.1, 0.2, 0.3, 0.4, and 0.5. 
(b) The peak-to-valley current ratio and (c) the biases for the 
current peak, kink, and valley points as a function of Ah. 
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FIG. 8. (color online) (a) I-V characteristics for the graphene 
DB RTD considered in Fig. 5 with 5 = 2.5 and different 
ADBWs. Along the arrow, A,,-0, O.f , 0.2, 0.3, 0.4, and 0.5. 
(b) The peak-to-valley current ratio and (c) the biases for the 
current peak, kink, and valley points as a function of A^,. 



calculations show that bigger PVRs can also be achieved 
by a joint of both asymmetries of heights and widths. 
However, the effect for a value of A^, -|- A^, is somewhat 
the same for a same value of an individual Xh or A^^. The 
different dependencies of PVR on band gap and ADBH 
or ADBW come from that the band gap not only sup- 
presses hole-to-electron transport and Klein paradox but 
also enhances resonant tunneling, while the structural 
asymmetries only enhance resonant tunneling in a proper 
range. 

In Sec. IV, we have seen that the first NDR operation 
window is pinked by the kink effect. Now Figs. 6(c)- 
8(c) tell us that, the kink effect is also robust against 
the band gaps or structural asymmetries. As a result, 
even with these factors, the first NDR operation window 
can be solely determined by the Fermi energy. This is 
a totally new feature for NDR in graphene that has not 
been found in any other types of graphene NDR devices. 



VI. CONCLUSIONS AND REMARKS 

In summary, we have constructed an analytic model 
which takes into account the realistic linear voltage drop 
and three transport regimes to investigate the NDR ef- 
fect of massless Dirac Fermions in graphene DB RTDs. 
The I-V and DC characteristics are calculated within 
the Landauer-Biittiker formalism in the rotated pseu- 
dospin space. The results indicate that, in graphene 
DB RTDs the NDR is a result of competition between 
hole-to-electron transport, Klein paradox, and resonant 
tunneling, and it only presents for appropriate structural 



the NDR strength by a tunable band gap in graphene 
and structural asymmetries of the DB RTD and investi- 
gated their roles on the PVR. Our results show that the 
PVR can be enhanced near exponentially by the induced 
band gap while it displays a parabolic-like dependence on 
the asymmetric degree in barrier heights and/or widths. 
We also observe a totally new and remarkable feature: 
the operation window of NDR in graphene DB RTDs 
is always pinned by a robust kink effect hence can be 
controlled easily by solely tuning the Fermi energy. This 
fact provides an opportunity for designing NDR based de- 
vices, such as high-speed bistable switching circuits and 
very high frequency oscillators, in a controllable way. 

Our calculations also show that the differential con- 
ductance displays totally different behaviors in the three 
transport regimes, and moreover, the NDR is absent in 
graphene SBDs which is in contrast to the earlier reports. 
Noting that the NDR only occurs in DB RTD with proper 
parameters, it is not surprising that the NDR never oc- 
curs in SBDs where no resonant tunneling takes place. 
Although we take graphene as an example, the results 
also apply to surface states of 3D topological insulators 
[3^ 1 where the real electron spin rather than the sublat- 
tice structure in graphene provides the Dirac structure. 

VII. ACKNOWLEDGEMENT 

This work is supported by the NSFC (10974109 and 
11174168), the SRFDP (20100002110079), and the 973 
Program of China (2011CB606405). HCW was grateful 
to the SFI Short Term Travel fellowship support during 
his stay at PKU. 



7 



[1] H. Mizuta and T. Tanoue, The Physics and Applica- 
tions of Resonant Tunneling Diodes. Cambridge Univer- 
sity Press, Cambridge, 1995. 

[2] W. Middleton, M.E. Van Valkenburg, Reference Data for 
Engineers: Radio, Electronics, Computers and Commu- 
nications. Boston Oxford, Newnes, 2002. 

[3] A. Abdel-All, A. Elshafie, and M.M. Elhawary, Vacuum 
59, 845 (2000). 

[4] J. McCinness, P. Corry, and P. Proctor, Science 183, 853 
(1974). 

[5] H. Shirakawa, E.J. Louis, A.G. MacDiarmid, C.K. Chi- 
ang, and A.J. Heeger, J. Chem. Soc, Chem. Commun. 
578 (1977). 

[6] D. Dragoman and M. Dragoman, Appl. Pliys. Lett. 90, 

143111 (2007). 
[7] V. Nam Do, Appl. Phys. Lett. 92, 216101 (2008). 
[8] V. Nam Do, V. Hung Nguyen, P. Dollfus, and A. Bournel, 

J. Appl. Phys. 104, 063708 (2008). 
[9] Z.F. Wang, Q. Li, Q. W. Shi, X. Wang, J. Yang, J.G. 
Hou, and J. Chen, Appl. Phys. Lett. 92, 133114 (2008). 
[10] H. Ren, Q.-X. Li, Y. Luo, and J. Yang, Appl. Phys. Lett. 

94, 173110 (2009). 
[11] H. Fang, R.-Z. Wang, S.-Y. Clien, M. Yan, X.-M. Song, 

and B. Wang, Appl. Phys. Lett. 98, 082108 (2011). 
[12] V. Nam Do and P. Dollfus, J. Appl. Phys. 107, 063705 

(2010) . 

[13] V. Hung Nguyen, A. Bournel, and P. Dollfus, J. Appl. 

Phys. 109, 093706 (2011). 
[14] V. Hung Nguyen, F. Mazzamuto, J. Saint-Martin, A. 

Bournel, and P. Dollfus, Appl. Phys. Lett. 99, 042105 

(2011) . 

[15] Gerson J. Ferreira, Michael N. Leuenberger, Daniel Loss, 

and J. Carlos Egues, Phys. Rev. B 84, 125453 (2011). 
[16] Y. Wu, D. B. Farmer, W. Zhu, S.-J. Han, CD. Dimi- 

trakopoulos, A. A. Bol, P. Avouris, and Y.-M. Lin, ACS 

Nano 6, 2610 (2012). 
[17] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, L V. Crigorieva, and A. A. 

Firsov, Science 306, 666 (2004). 
[18] T.C.L.G. Sollner, W.D. Goodhue, P.E. Tannenwald, 

CD. Parker, and D.D. Peck, Appl. Phys. Lett. 43, 588 

(1983); A.R. Bonnefoi, R.T. Colhns, T.C. McGill, R.D. 

Burnham. and F.A. Ponce, ihid. 46, 285 (1985). 
[19] B. Huard, J. A. Sulpizio, N. Slander, K. Todd, B. Yang, 

and D. Goldhaber-Gordon, Phys. Rev. Lett. 98, 236803 

(2007); S. Russo, M. F. Craciun, M. Yamamoto, S. 

Tarucha, and A. F. Morpurgo, New J. Phys. 11, 095018 

(2009). 

[20] J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and 
C.W.J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 

[21] The kink effect is an interesting phenomena that the min- 
imum differential conductance happens exactly at a bias 
equals to the Fermi energy and the operation window 
is the bias range between the current peak and valley 
points, where DC keeps negative. 

[22] A. Javey, H. Kim, M. Brink, Q. Wang, A. Ural, J. Guo, 
P. Mcintyre, P. Mceuen, M. Lundstrom, and H.J. Dai, 
Nature Mater. 1, 241 (2002). 

[23] Y. Song, H.C Wu, and Y. Guo, Appl. Phys. Lett. 100, 
253116 (2012). 

[24] Y.M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 



[25] Y.-B. Tang, et al. ACS Nano 7, 1970 (2012). 

[26] Y. Song and Y. Guo, J. Appl. Phys. 109, 104306 (2011). 

[27] M. Born and E. Wolf, Principles of optics: electromag- 
netic theory of propagation, interference and diffraction 
of light. Oxford, Pergamon Press, 1964. 

[28] I. Meric, M.Y. Han, A.F. Young, B. Ozyilmaz, P. Kim, 
and K.L. Shepard, Nature Nanotech 3, 654 (2008). 

[29] D. H. Chow and T. C. McGill, Appl. Phys. Lett. 48, 1485 
(1986); D. H. Chow, T. C McGill, I. K. Sou, J. P. Faurie, 
and C. W. Nieh, ibid. 52, 54 (1988). 

[30] R. Beresford, L. F. Luo, and W. I. Wang, Appl. Phys. 
Lett. 54, 1899 (1989); J. R. Soderstrom, D. H. Chow, 
and T. C. McGill, ibid. 55, 1348 (1989). 

[31] P.N. Racec, T. Stoica, C. Popescu, M. Lepsa, and Th. 
G. van de Roer, Phys. Rev. B 56, 3595 (1997). Since the 
right barrier will be lower more than the left barrier under 
the bias, a right barrier with larger height or bigger width 
(i.e., positive or A^) should be excepted to obtain 
stronger NDR. 

[32] X.L. Qi and S.C Zhang, Rev. Mod. Phys. 83, 1057 
(2011). 



